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Abstract 

We study the effect of an externally imposed oscillatory shear on the motion of a grain bound- 
ary that separates differently oriented domains of the lamellar phase of a diblock copolymer. A 
direct numerical solution of the Swift-Hohenberg equation in shear flow is used for the case of a 
transverse/parallel grain boundary in the limits of weak nonlinearity and low shear frequency. We 
focus on the region of parameters in which both transverse and parallel lamellae are linearly stable. 
Shearing leads to excess free energy in the transverse region relative to the parallel region, which 
is in turn dissipated by net motion of the boundary toward the transverse region. The observed 
boundary motion is a combination of rigid advection by the flow and order parameter diffusion. 
The latter includes break up and reconnection of lamellae, as well as a weak Eckhaus instability 
in the boundary region for sufficiently large strain amplitude that leads to slow wavenumber read- 
justment. The net average velocity is seen to increase with frequency and strain amplitude, and 
can be obtained by a multiple scale expansion of the governing equations. 
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I. INTRODUCTION 



Below an order-disorder temperature, nanoscale phases of various symmetries can be 
found in block copolymer melts For example, for diblock copolymers that consist 

of two chemically distinct but covalently bonded monomers, six distinct phases that result 
from microphase separation have been documented jll . They differ by the symmetry of their 
composition modulation. One of the most actively investigated microphases is the lamellar 
phase, which can be observed around 50-50 composition (symmetric mixture). However, 
when melts are processed by thermal quenching or solution-casting from a disordered phase, 
a macroscopic size sample usually exhibits polycrystalline configurations comprised of many 
locally ordered but randomly oriented domains (grains) and large amounts of topological 
defects, such as grain boundaries, dislocations, and disclinations. Such state of partial 
ordering is undesirable for many applications, with significant affects on e.g., the optical 
and mechanical performance of the material. Therefore, the realization and control of long- 
range orientational order of domains has been of great interest in both experimental and 
theoretical studies. 

Several methods have been used experimentally to achieve microstructural alignment in 
block copolymer melts since the discovery of flow- induced alignment by Keller et al. . One 
of the most common ways to induce global order in bulk samples is to impose an oscillatory 
or steady shear between two parallel plates, which not only has the advantage of easily char- 
acterizing the shear stress and monitoringthe aligning progress, but also can bring some 
interesting new physics . For the case of lamellar morphology that we 

are interested in here, the two most often observed alignments relative to the imposed shear 
are parallel (with the lamellar layers parallel to the shearing plane), and perpendicular (with 
the lamellar layers normal to the vorticity of the shear flow). Both alignments have been 
observed in nearly symmetric diblock copolymer systems such as poly(ethylenepropylene)- 
poly(ethylethylene) (PEP-PEE) and polystyrene-polyisoprene (PS-PI) 0, Which of 
these two possible alignments is selected in a given experiment depends on temperature, and 
shear strain amplitude and frequency ^ 0, |^ . A third alignment direction, the so called 
transverse orientation in which the lamellar normal is parallel to the shear direction, has been 
found to coexist with parallel orientation for entangled poly(styrene-b-ethylenepropylene) 
(S-EP) diblock copolymer at high frequency and in the strong segregation limit IQj. Other 
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studies have focused on the kinetics of global alignment and have shown that the ordering 
rate increases with shear frequency, strain amplitude, and temperature Experiments 
on PS-PI diblocks further show nonlinear effects of the strain amplitude on the alignment 
rate, and that the time scale for the development of alignment exhibits a power law de- 
pendence on the strain amplitude, with an exponent equal to —3 or —5 depending on the 
stage of alignment and in different frequency regimes ll|. Furthermore, many experiments 
have indicated that the motion of topological defects plays an important role in the global 
alignment of microdomains under shear, including the evolution of kink band defects and 
tilt boundaries [3], as well as the migration and annihilation of partial focal conic defects, 



boundaries, and tilt walls 



A coarse-grained description of a diblock copolymer melt has been used, both in the weak 
segregation 
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115|,Q egation [15[ regimes, to theoretic ally understand the 

effect of shear flow and the mechanisms of morphology evolution 

liaHHQ- The 

analytic studies of Gates and Milner ^3], and Fredrickson focused on the order- disorder 
transition from the isotropic state to the lamellar phase and the related alignment dynamics 
for systems subjected to steady shear flow. Below the transition temperature, the stability of 
uniform lamellar structures under an oscillatory shear flow has been addressed in both two- 
dimensional (2D) and three-dimensional (3D) [2flj cases. Bifurcation diagrams, including 
secondary instabilities (Eckhaus and zigzag) were given. These analyses found that although 
all three types of lamellar orientations could be linearly stable under specific conditions, the 
stability range of the perpendicular orientation is larger than that of the parallel one, which 
in turn is larger than the region of stability of the transverse orientation. The perpendicular 
alignment was also shown to be the preferred orientation following the decay of unstable 
parallel or transverse lamellae. In addition to these analytic work, computer simulations 
have been performed to investi gate the dynamics of lamellar alignment in bulk samples 
under steady or oscillatory [2^ shear flow. The effects of strain amplitude and shear 
frequency on the degree of sample alignment have been examined, as well as annihilation 
processes of defects (such as dislocations and disclinations) [2^. Also, domain coarsening in 
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2D diblock copolymer melts has been addressed although in the absence of shear 
Little is known about the effect of shear on coarsening of block copolymers. 

Few of the previous theoretical studies discussed above focused on the detailed dynamics 
and quantitative properties of topological defects motion under shear, which are crucial for 
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the understanding of alignment and coarsening. In this paper we study the detailed mech- 
anisms of grain boundary motion under an oscillatory shear flow building upon stability 



results of uniform lamellar patterns 
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Is^. We use a simplified 2D configuration which 
involves only two lamellar domains of parallel and transverse orientations respectively, sep- 
arated by a grain boundary, and focus on shears of small amplitude (less than 50%) and 
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low angular frequency. Compared to previous numerical simulations 
aspect ratio (defined as the ratio between the lateral extent of the system and the lamellar 
wavelength) is larger, so that important dynamic features associated with grain bound- 
ary motion, such as diffusive relaxation of lamellae, phase shift of boundary velocity, and 



wavenumber adjustment in the transverse region, can be quantitatively analyzec . 
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features were absent in earlier work on lamellar (roll) systems without shear 
and we argue below that they are important when the system is under oscillatory shear. 

This paper is organized as follows: In Sec. ^Iwe introduce a dimensionless model equa- 
tion based on the Swift-Hohenberg equation to describe the dynamic evolution of symmetric 
diblock copolymer melts under shear, and describe the grain boundary configuration used. 
The numerical results including grain boundary velocity and lamellar wavenumber are pre- 
sented in Sec. IIIII We derive amplitude equations governing the system evolution in Sec. 
rvl and compare the results with the direct solutions of Sec. IIIII Finally, in Sec. El we 
summarize our results and discuss the physical origin of the phenomena observed. 



II. MODEL EQUATION AND GRAIN BOUNDARY CONFIGURATION 

The system under consideration is a symmetric diblock copolymer melt below the 
order-disorder transition temperature Tqdt- For length scales larger than the microscopic 
monomer scale (i.e., at a mesoscopic level) and time scales long enough compared to the 
molecular relaxation of the polymer chains, a coarse-grained description of the block copoly- 
mer melt can be used, with an order parameter field ip{r,t) representing the local density 
difference of the two constituent monomers. In the weak segregation limit, i.e., close to 
Tqbt, a coarse-grained free energy functional has been derived 14, 15, Q: 
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where r denotes a reduced temperature variable which is a measure of the distance from the 
order-disorder transition and is positive below Tqdt, and = 2ti /X^ is the wavenumber of 
the periodic lamellar structure. Under the assumption that changes in the local composition 
field if) are driven by the free energy minimization and advection by the flow, the order 
parameter '?/'(r, t) obeys a time-dependent Ginzburg-Landau equation 

^+,.V^ = -A-. (2) 

Here v is the local velocity field, A is an Onsager kinetic coefficient and can be written as 
A = Mqf (with M a mobility). 

We introduce a length scale I/qq, a time scale I/A^Qq^ (which is the characteristic polymer 
relaxation time and ~ l/Dq^'^, with D the chain diffusivity of the copolymer), and an order 
parameter scale (^q^'^/uY^'^. Given PEP-PEE-2 as an example, we have Aq ~ 30 nm, 
D ~ 10^^^ cm^/s for temperatures close to Tqdt; and hence the length scale here is about 5 
nm and the time scale is about 0.03 s. Consequently, the rescaled composition field obeys 
a dimensionless Swift-Hohenberg equation j30| with an advection term 

^ + V ■ Vz^ = ez^ - (V^ + go') V - (3) 

where e = t/^^q^ and we have < e ^ 1 in the weak segregation regime considered here. 
Also, qo = 1 although the symbol go is retained in what follows for clarity of presentation. 

As shown in Fig. ^ we consider a 2D reference state below the order-disorder transition 
containing a planar grain boundary that separates two semi-infinite ordered domains A and 
B. Initially both domains are in the lamellar state with the same wavenumber go but oriented 
along different directions. We are interested here in the case of a 90° grain boundary with 
two mutually perpendicular lamellar sets A and B, a co nfigu ration that is known to be stable 



against small perturbations in the absence of shear 
an imposed shear fiow 
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The two domains are under 



vo = "^^^ = 7^ cos(ti;t) z X, (4) 

where da/ dt represents the shear rate with strain a{t) = 7 sm{ujt), uo is the angular frequency, 
7 the strain amplitude, and all quantities are assumed dimensionless. Lamellae of domain 
A are transverse (with the wavevector components q^ = go and g^ = 0) at t = 0, and those 
in region B parallel {qx = and qz = go)- Parallel lamellae B are marginal to the shear 
and not distorted, while transverse lamellae A are compressed, with both orientation and 



wavelength changing following the imposed shear as shown schematically in Fig. Thus, 
we anticipate that the grain boundary will not remain stationary even though both A and B 
lamellae are linearly stable under shear. Net motion results from the free energy difference 
between region A (compressed) and region B (unchanged), as well as diffusive relaxation of 
the order parameter as shown below. 

The stability of a uniform co nfig uration of either parallel or transverse lamellae under 



shear flow has been given in Refs. 



19 and 



20l . There exists a critical strain amplitude 7c above 



which the lamellar structure of a given orientation melts, with 7^ — * oo for parallel lamellae 
of wavenumber q = qo, and small 7c for transverse orientation. The stability diagrams 
presenting secondary instability boundaries (for zigzag and Eckhaus modes) for 2D system 



have also been given in Ref. 



I9I . We focus below solely on shears for which both uniform 



parallel and transverse lamellae are linearly stable. In addition, we consider the case in which 
shear effects are of the same order of magnitude as diffusive relaxation of the order parameter. 
Otherwise, at one extreme lamellae are passively advected by the flow, whereas at the other, 
diffusion dominates. If the velocity v in Eq. ^ can be approximated as vq (which is the 
case under certain conditions including neglecting back flows due to osmotic stresses, 
and any viscosity contrast between the microphases) , then the interesting range of 7 is such 
that the advection contribution due to the imposed shear (v ■ V = {da/dt)zdx in Eq. Q) is 
0{e). As will be further discussed in Sec. IIVI in connection with our multiple scale analysis, 
this requires 7 ~ 0{e^^^), an assumption that will be used in what follows. 



III. NUMERICAL RESULTS 

We first introduce a time dependent sheared frame of reference 3, in which the 
imposed shear flow vanishes. It is defined by the non-orthogonal basis set {e^i = x, e^/ = 
a(t)x + z} shown in Fig. ^ In this sheared frame, we have the coordinates x' = x — a(t)z and 
z' = z. Also, the corresponding reciprocal basis set is defined as {g^' = x — a(t)z, g^/ = z}, 
with wavevector components expressed as q^' = qx and g^/ = a(t)qx + qz- After coordinate 
transformation the Swift-Hohenberg equation Q becomes 

^=e^- (V^ + qin - ^^ (5) 
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where the modified Laplacian operator takes the form 

V'^ = [1 + a\t)] dl - 2a{t)d^>d,^ + dl. 

The critical value 7c for neutral stability is independent of the shear frequency u and given 
by 



-b + Vb^ - 4:dc 



^ ) ' 

with b = ql{2(3^ + Aql - 2ql)/2, c = (3^- 2go/?' + - e, and d = 3g^,/8(here = ql + ql). 
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The results for secondary instabilities were already presented in Ref. 
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We numerically solve Eq. (0) with periodic boundary conditions (in the sheared frame) by 
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and also detailed in the appendix of 



usin g th e pseudo-spectral algorithm described in Ref. 
Ref. llSl. The equation has been discretized on a 1024 x 1024 square grid in the sheared frame. 
In most of our calculations, we choose a mesh spacing Ax' = Az' = Ao/8, corresponding to 
8 grid points per wavelength Aq (= 27r/go) and so Aq = 1/128 for the wavenumber spacing. 
A second order, semi-implicit time stepping algorithm is used, with a time step At = 0.1. 
In order to remain within the weak segregation regime we set e = 0.04 and the angular 
frequency u chosen is 0{e). Initial conditions are obtained by numerical solution of Eq. (0) 
without shear, i.e., a{t) = 0, from a starting configuration consisting of two symmetric grain 
boundaries located at x' = Lx'/4: and 3Lx'/4: (with L^' = 1024 Ax' the extent of the system in 
the x' direction), separating a parallel domain B from two surrounding regions of transverse 
lamellae. This configuration is allowed to evolve without shear until a stationary solution 
is reached. This stationary solution is used as the initial condition for the integration with 
a{t) 0. 



A. Grain boundary motion due to shear and diffusive relaxation of the order 
parameter 

If 7 > 7c = (8e/3)^/^, the stability limit of transverse lamellae (Eq. ©), A lamellae melt 
to a disordered state {ip = 0), and B lamellae invade region A. If 7 < 7c, but still large 
enough to result in a secondary instability (Eckhaus, cross-roll, or zigzag) of transverse 
lamellae at a given frequency (e.g., 7 = 0.5 for e = 0.04 and small u), our calculations 
show that small domains of parallel lamellae form within the bulk transverse region A as a 
result of the instability. These domains evolve, and connect with each other and with the 
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approaching grain boundary to increase the extent of the parallel region B. In both cases, 
the final configuration of the system is a uniform lamellar structure of parallel orientation. 
More interesting phenomena are observed in the range of 7 and uo in which both parallel and 
transverse bulk regions are linearly stable, and when the contributions from shear flow and 
order parameter diffusion are of the same order, as described at the end of Sec. EI A typical 
transient configuration in this parameter range is shown in Fig. |21 with 7 = 0.4, uj = 0.04, 
and grid spacing Aa; = Ao/32. The configuration shown corresponds to t = 2.25To, where 
Tq = 2tt/uj is the shear period, and is presented in the laboratory frame basis set {x, z}. 

Our analysis of the transient evolution of the configuration under shear is based on the 
average location of the grain boundary x'^^. To determine this quantity we use a relation 
similar to that of Ref . hd B{x'^^) = ^3 EILJV^Kb, -^Kb. (i - l/2)Ao)]/4n = 6, with 
n the number of pairs of lamellae in the z' direction, and 6 a quantity 0{e). The value used 
here is 6 = e/4. The grain boundary velocity v'^-^ is defined as the time rate of change of 
Xgt,. Representative results (in the sheared reference frame) for x'^^ and function 
of time are shown in Figs. El and El Two distinct features can be clearly distinguished as 
illustrated in Fig. |21 First, transverse lamellae in the bulk rigidly follow the oscillatory shear 
fiow, leading to a periodic change in their orientation. Second, during part of the cycle the 
region of parallel lamellae B grows as transverse lamellae in the boundary region break up 
and reconnect as parallel lamellae (cross roll instability). This process is partially reversed 
during the rest of the cycle. The grain boundary exhibits oscillatory motion with a nonzero 
net average as shown in Figs. El (location) and 01 (velocity). Those portions of the shear 
cycle in which broken transverse lamellae recombine correspond to the segments in Fig. El 
with decreasing Xgj,, or negative velocity fgj, in Fig. jH 

Whereas rigid distortion of transverse lamellae is the dominant response to shear in the 
bulk, order parameter diffusion is important in the boundary region, and is the mechanism 
that enables dissipation of the excess free energy stored in the bulk transverse region due 
to shearing. On the one hand, transverse lamellae are elastically compressed by the shear, 
resulting in a net free energy increase in region A relative to B. As a consequence, the elastic 
contribution to the system's energy is expected to drive grain boundary motion toward the 
transverse domain at all times during the shear cycle. On the other hand, backward motion 
is also observed in Figs. ElandjUin those portions of the shear cycle in which the magnitude 
of the shear strain a{t) = 7sin(u;t) is small, indicating diffusive relaxation of the order 
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parameter field near the grain boundary. The competition between the two determines the 
net rate of advance of the boundary. This competition is illustrated in Fig. |^b). For 7 = 0.3 
and 0.4 (with the same frequency uo = 0.01), the figure shows that within one period Tq, a 
smaller strain amplitude 7 corresponds to larger temporal range of negative velocity, as well 
as to a larger phase lag between the boundary velocity and the imposed shear. The same 
effect can be seen in Fig. EJa), as lowering the frequency enhances the effect of diffusion 
over elasticity. 

The net velocity of the average location of the boundary is positive as shown in Figs. El 
andini The figures plot the temporal average of velocity over a period (fgj^) = dtv'^^{t) /Tq 
as a function of uj and 7. The velocity increases sharply at very small uj and saturates at 
large uj. Diffusive relaxation is more pronounced at lower u;, consistent with the calculations 
shown in Fig. Ufa), {v'^,^) also depends on strain amplitude 7, as seen in Figs. El and 
The average boundary velocity increases approximately as (fgb) ~ 7" with a ~ 4. The 
range of strain amplitude accessible to our calculations is limited by the restriction that 
7 ~ 0{e^^^). As discussed above, the transverse lamellar region is unstable for larger 7, 



N 



ote that the power law 



whereas for smaller 7, diffusion of order parameter dominates. 

dependence is consistent with experiments in PS-PI copolymers js, IjJ, in which the rate of 
global alignment of a bulk sample is a power law of the strain amplitude, with an exponent 
in the range 3 — 5. 



B. Wavenumber adjustment in the transverse region 

Bulk transverse lamellae elastically compressed by the shear would have a wavenumber 
go\/l + a?{t) in the laboratory frame, or q^i = qo, constant, in the sheared frame. Order 
parameter diffusion in the boundary region, however, is seen to lead to a wavenumber mod- 
ification of transverse lamellae relative to what would be expected from rigid deformation. 
We show in Fig. [Tithe structure factor l^'g^/l along the x' direction defined as the Fourier 
transform of the order parameter %p at fixed z' (e.g., at z' = Lz'/2, with L^' = 1024A2;' the 
system size along z' direction). For large enough 7 (7 > 0.1), the peak in lipg^,] shifts away 
from go = 1 with time, asymptotically reaching a constant < go (Fig- (Zfa)). We have 
studied this wavenumber compression for different frequencies u and strain amplitudes 7. 
We find that the value of is independent of uj, but that it increases with decreasing 7 
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(Fig. Htb)): = - = 5Aq (with Aq = 1/128) for 7 = 0.4 (solid line), 5q^ = 3Aq 
for 7 = 0.3 (dotted line), 6q^ = Aq for 7 = 0.2 (dashed line), and Sq"^ = for 7 = 0.1 (thin 
solid line). These values correspond to the disappearance of 5, 3, 1, and lamellae in region 
A respectively. We further discuss this finding in Sec. IIVI and in the discussion section IVl 
In addition to the motion in the position of the peak of the structure factor. Fig. [7|^a) also 
shows a decreasing amplitude and broadening of the peak. This is due to the finite size of 
the configuration; as the grain boundary moves, the region occupied by transverse lamellae 
decreases. 

The power spectrum of region B {ipg^, \ as a function of qz' (the Fourier transform of ip at 
fixed x' position, results not shown here) is unaffected by the shear flow. Its maximum is 
located at q^ = q^i = qo, without any visible shift in time. 



IV. MULTIPLE SCALE ANALYSIS AND AMPLITUDE EQUATIONS 



A multiple scale analysis of the 
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type frequently used to derive amplitude equations close 



to instability thresholds j26 



32l | can be introduced here to derive an equation of motion 
for the grain boundary in the limit of weak segregation e <^ 1. As shown in Fig. ^ we 
flrst deflne a time dependent, orthogonal basis set {e^^ = (x — a(t)z)/{l + a^(t)), e^^ = 
(a(t)x + z)/(l + a^(t))} (different from the non orthogonal sheared frame {e^^', e^} used 
in Sec. II11|) . This frame is the orthogonal frame of reference attached to the transverse 
region. The corresponding coordinates are xa = x — a(t)z and = a(t)x + z. We then 
introduce an anisotropic coordinate scaling to deflne slowly varying amplitudes of exp{iqoXA) 
as Xa = e^^'^XA and Za = e^^^ZA- We retain the laboratory frame coordinates in region B 
{x, z} with a base mode given by exp{iqoz). Its slowly varying amplitude is a function of 
the rescaled variables Xb = e^^'^x and Zb = e^^'^z. We then expand the order parameter 
fleld ip as 

= (Ae^""^^ + Be'^^' + c.c.) , (7) 

where both complex amplitudes A and B are functions of the slow spatial scales Xa, Za, 
Xb, and Zb, and of a slow time scale T = et. The requirement of slow amplitude change 
restricts our analysis to low frequencies u ~ 0{e), and we further focus on sufficiently 
small shear amplitudes so that advection and local diffusion balance. This requires v ■ V = 
{da/dt)zdx ~ (V^ + 90)^ ~ ^ according to Eq. (jHl). Considering that in the transverse region 
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dx = dxA + = -<^dxA +dzAi z = {-axA + -2a)/(1 + a^), as well as the spatial {xa, za) 

and temporal scalings, we require that 7 ~ 0{e^^'^). The same relationship follows from the 
scahngs appropriate for the parallel region. 

Following standard multiple scale procedure ) we introduce the expansions — *■ 



XA 



'/\dxA + adzA), d,^d, + e^/^-adxA + ^^J + e^^dz, - e^/^a^x,, and 
dt + {da/dt)zdx t[dT + '^n{XAdzA + ZBdxg)] + 0{e'/^) (here d and 7^ are defined by 
a = e^^^d and da/dt = e^^^jn), and we derive the following amplitude equations at (9(e^/^) 
from Eq. Q, 



dtA 
dtB 



da 



e - ^XAdzA - i'^iqodxA + " ^o^^Y 



dt 
da 
'dt 



A - UrA-2|5rA 



B- \B\'^B -2\A\'^B. 



(8) 
(9) 



28. 



For 7 = 0, these equations reduce to those of Refs. |27| and 

Equations (jH]) and (jH]) are expressed in two different coordinates systems. We next trans- 
form them to a common sheared frame {e^/, e^/}, by using the relations xa = x' and 
Za = CLx' + (1 + a'^)z'. The base state of the order parameter is still given by Eq. (|7j) (with 
Xa replaced by x'), and to (9(e^/^), the resulting 2D amplitude equations were already given 



in Ref. 



33I . Here we further assume a planar grain boundary in the sheared frame, for which 



the dependence of the amplitudes on the coordinate {z') parallel to the grain boundary can 
be ignored. The complex amplitudes A and B satisfy the one-dimensional (ID) equations 



dtA =[e- {2iqodx' - qla^] A - \A\^A - 2\B\^A, 



(10) 



and 



dtB =[e- {-2iqoadx' + d^,)^] B - \B\'^B - 2\A\'^B. (11) 

At (9(e^/^), two contributions from the shear flow remain in Eqs. ()10|) and pij) . The first 
one involves the term —2iqQadx' in Eq. (fTTj) and the term (gQa^)(2igo<9x') obtained from the 
expansion of Eq. (fTUI) . which is non negligible only in the grain boundary region and leads 
to diffusive relaxation of the order parameter. The second is q^a^A in Eq. (jiop . This term 
is uniform in the entire region A, and reflects the contribution from advection of transverse 
lamellae by the flow. 

We had analytically calculated the velocity of the grain boundary from these amplitude 
equations in Ref. by assuming that the amplitudes can be approximated by A{x',t) ~ 
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A{x' — x'g^(t)) and B{x',t) ~ B{x' — Xgj,(t)). We found that the velocity is proportional 



to the free energy difference between t 
previous studies in the absence of flow 



le transverse and parallel phases, in agreement with 



2i 
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29|. Also, the results gave the correct order of 



magnitude of the average velocity, but we noted quantitative discrepancies j33|. We argued 
that the adiabatic approximation for A{x',t) and B{x',t) given cannot incorporate diffusive 
relaxation of the order parameter in the boundary region so that the calculation only yields 
an upper bound to the net boundary velocity. Since we have argued in Sec. IIIII that this 
relaxation is important, we turn here to a numerical determination of the boundary velocity. 

We present next the results of the numerical solution of Eqs. and (fTT|) . We initially 
consider a region of parallel lamellae B surrounded by two regions of transverse lamellae 
A, and use periodic boundary conditions in the integration. Both amplitudes A and B are 
complex variables. A pseudo-spectral method is applied, with a Crank-Nicholson scheme 
used for the linear terms, and a second-order Adams-Bashford scheme used for the nonlinear 
terms. The instantaneous location of the grain boundary x'g-^^{t) is defined by the condition 
l-^(^gb)l — and its velocity v'^^ is found by taking the time derivative of x'^^{t). In 
order to compare with the 2D results of the original model shown in Sec. IIIH we set the 
system size L = 1024, the time step At = 0.1, and the grid spacing Ax' = Ao/8. As 
was the case there, the initial condition for A and B is provided by the steady solution of 
the amplitude equations in the absence of shear. Our results for the boundary velocity for 
7 = 0.4, to = 0.01, and u = 0.04 are shown as dotted curves in Fig. Ufa), both in good 
agreement with the direct solution of original model equation (0) (symbols in the figure). 
The time averaged velocity (fgi,) is shown by the dashed lines in Figs. ElandlHl 

In order to further analyze the wavenumber readjustment process discussed in Sec. IIII B| 
we show our results in terms of the phase 0^ of the complex amplitude A. In the sheared 
frame we define A = \A\ exp(i0^) and plot (pA as a function of grid index i^' = x' /Ax' in 
Fig. IHl Near the boundary, the phase becomes linear in space: (pA oc —6q ■ x' (the linear 
behavior is clearer for a larger system size, as seen by comparing Fig. IHfa) (L = 1024) 
withlHfb) (L = 4096)), indicating a local wavenumber change Qx' — > Q'o — Sq- This is also in 
agreement with the direct solution of the original model as shown in Fig. [7| Note also that 
the region of linearity (right side of dot-dashed line in Fig. |H)) increases with time, indicating 
that the re-adjustment of the local wavelength of the transverse lamellae first occurs at the 
boundary, and then progressively propagates into the bulk. At late times (e.g., t = 40To in 
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Fig. IHl^a)), the wavenumber change can be observed in the whole transverse domain, with 
g^/ corresponding to the stationary peak position of the structure factor \ipq^i \ presented in 
Fig. [7|^a) {t > 40To there). We find that 5q = 5q^, the wavenumber shift discussed in Sec. 
11111 by determining 5q from the slope of the dotted line in Figs. |Hl^a)and (b). 

V. DISCUSSION AND CONCLUSIONS 

A coarse grained order parameter model has been used to study the motion of a grain 
boundary separating two regions of uniform parallel and transverse lamellae under an im- 
posed shear flow. The motion of the boundary is oscillatory, and the driving force for motion 
is the excess energy stored in the elastically strained transverse phase that can only be re- 
lieved through diffusive relaxation of the order parameter in the boundary region. Diffusive 
relaxation, however, is complex as the response of the order parameter field is out of phase 
with the shear, and lamellae break up and reconnect during each of the cycles. As expected, 
the effects of diffusive relaxation are more pronounced for small shear strain and low angular 
frequency, as seen in both the time dependent behavior of boundary velocity (Fig. or 
the time averaged velocity (Figs. ElandlHl). Although under the conditions of the study 
both transverse and parallel orientations are linearly stable, we observe net motion of the 
boundary toward the region occupied by transverse lamellae. 

As the boundary moves over time, we have observed that the wavenumber of the trans- 
verse lamellae does not remain constant and equal to go- Instead, it is slowly readjusted 
by wavenumber diffusion, as shown by both direct solution of the coarse grained model of 
Sec. mil and the corresponding complex amplitude equations of Sec. |TVl The wavenumber 
shift 5q^ is approximately independent of shear frequency cj, but strongly dependent on the 
strain amplitude 7. In order to understand the physical origin of wavenumber compression 
(or lamellae expansion) that occurs in the transverse region A, we focus on the lamellae 
near the grain boundary, since the calculation in Sec. II VI shows that wavenumber change is 
initiated at the boundary, and then propagates into the bulk. Since the amplitude of the 
transverse lamellae go to zero at the boundary region, it is possible to create or eliminate 
lamellar planes there in a way that it is not possible in the bulk for the parameters of our 
study I19I . First consider the stability diagram of the Swift-Hohenberg model at zero 
shear jSJ]. For fixed e and q^ > qo, the closest instability boundary in the q^-e diagram is 
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the Eckhaus instability given by 

e = 12(g, - go)'. (12) 

For e = 0.04 and go = 1, we have the Eckhaus boundary at g;^ = go + (e/12)^/^ = 1.0577. 
In our case g^ is the wavenumber of the transverse lamellae in the lab frame, and equal 
to qo^/l + a'^{t) if we assume that the lamellae are rigidly distorted by the shear. The 
maximum value of qx is then g™^'' = go a/1 + 7^ so that for 7 = 0.4 we have q^^"^ > q^. 
Although this is not sufficient to destabilize the bulk transverse region (according to the 
stability diagrams in Ref. 1191 obtained by Floquet analysis over the entire period of the 
oscillation), it appears to be sufficient to induce lamellar elimination at the grain boundary 
region, as seen in Figs, d and |H1 The figures show that 5g^ = 5Ag, corresponding to the 
elimination of 5 transverse lamellae. As 7 decreases g™*^^ also decreases, becoming smaller 
than q^, and eventually lamellae elimination is expected to cease. This is consistent with 
our results shown in Fig. [Z(b). With decreasing value of 7 from 0.4 to 0.1, the number of 
lost transverse lamellae decreases from 5 to 0. 

In summary, for small shear strains and low frequencies such that diffusion of order pa- 
rameter is of the same order as advection by the flow, the excess free energy in the transverse 
region relative to the marginal parallel region is dissipated through order parameter diffusion 
in the grain boundary region. The latter includes break up and reconnection of transverse 
lamellae, and a weak Eckhaus instability developing at the grain boundary that diffuses 
into the bulk transverse lamellae leading to dynamical wavenumber readjustment. A weakly 
nonlinear analysis, as well as the amplitude equations derived, capture quantitatively all 
aspects of grain boundary motion, including the boundary velocity and the wavenumber 
readjustment. The order parameter distribution in the boundary region can be represented 
crudely by introducing an adiabatic approximation into the amplitude equations, which gives 
a reasonable approximation to the net boundary velocity toward the transverse region, and 
be well reproduced by the direct solution of the amplitude equations. Although our study 
is conflned to the case of a transverse/parallel grain boundary in two dimensions, we expect 
that our results will qualitatively hold for both three dimensional transverse/parallel and 
transverse/perpendicular cases. In three dimensions, however, there is a completely different 
type of tilt boundary, that between parallel and perpendicular lamellae. Both orientations 
are marginal with respect to the shear. This configuration is currently under investigation. 
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FIG. 1: Schematic representation of the two dimensional grain boundary configuration under 
oscillatory shear flow studied in this paper. Both the non orthogonal sheared frame {e^'je^'}, and 
the auxiliary frame {e^^^je^^} are indicated. 
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FIG. 2: Grain boundary configuration (in gray scale) at time t = 2.25To obtained by numerically 
solving the Swift-Hohenberg equation Q with e = 0.04, 7 = 0.4, uj = 0.04, and Ax = Ao/32. 
The configuration shown here has been transformed back to the laboratory frame {x, z}. As time 
evolves, the transverse domain is invaded by parallel lamellae, until a uniform parallel configuration 
occupies the whole system. 
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FIG. 3: Relative grain boundary displacement in the sheared frame {x'g^ — x'g^(t = 0))/Ao as a 
function of time t/To for e = 0.04, 7 = 0.4, u = 0.04, and Ax = Ao/8. 
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FIG. 4: Boundary velocity v'^^ as a function of time obtained from direct numerical integration 
of Eq. ((S)) with e = 0.04 and different values of uj and 7. (a) Fixed strain amplitude 7 = 0.4, 
iv = 0.01 (filled triangles), and uj = 0.04 (open circles). The dotted lines are the corresponding 
results obtained from numerical integration of the amplitude equations (fTU]) and (fTT|) . (b) Results 
at constant frequency to = 0.01, but different strain amplitudes 7 = 0.4 (solid line), and 0.3 (dashed 
line) . 
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FIG. 5: Temporal average of the grain boundary velocity {v'g^) as a function of frequency uj. The 
symbols correspond to the solution of the order parameter model Eq. (jSJ for 7 = 0.4 (circles), and 
0.3 (squares), while the corresponding dashed lines are obtained from the amplitude equations (|10jl 
and (HU). 
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FIG. 6: Log- log plot of average velocity ('Vgb) versus strain amplitude 7 for = 0.01 (squares), 
0.02 (circles), and 0.1 (triangles) from direct solution of Eq. ©. Also shown (dashed lines) are 
the corresponding results given by the amplitude equations (|lfl|) and (|1H) . 
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FIG. 7: One dimensional structure factor for transverse lamellae along the x' direction at 

z' = L^i jl as a function of wavcnumbcr q^i ■ (a), different times from top to bottom: t = 30To, 
40ro, SOTo and 60To, with 7 = 0.4 and u = 0.04. (b), different strain amplitudes: 7 = 0.4 and 
t = lOOTo (solid curve), 7 = 0.3 and t = SOOTq (dotted curve), 7 = 0.2 and t = 1200ro (dashed 
curve), and 7 = 0.1 and t = 12000ro (thin solid curve). Here u = 0.1, and At = 0.1 except for 
7 = 0.1 in which case we have used At = 0.2. In both cases (a) and (b), the vertical dot-dashed 
line indicates the location of the wavenumber q^' = qo = 1. 
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FIG. 8: Phase (pA of the complex amphtude ^ as a function of position in the sheared frame (in 
terms of the grid index ij./) with Ax' = Ao/8, 7 = 0.4, and uj = 0.04. Two system sizes are shown: 
(a) L = 1024 at times (from left to right): t = 2Tq, 5Tq, IOTq, 20Tq, and 40To (to be compared with 
Fig. Eta)); and (b) L = 4096 at times (from left to right) t = WTq, SOTq, IOOTq, and 200ro. The 
vertical dot-dashed lines indicate the instantaneous grain boundary positions x'^y^, and a dotted 
line with slope —SqZi ■ Ax' = — 5Ax'/128 is also shown for reference. 
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